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Velocity and two-directional focusing of charged particles 
in crossed electric and magnetic fields. I 


By Nits SvARTHOLM 


1. Introduction 


In all systems intended for focusing of charged particles and giving energy 
or momentum dispersion, one has always to compromise between the demand 
of high resolving power and that of a high transmitting power. When one 
goes beyond the most elementary focusing systems, e.g. those achieved by 
homogeneous magnetic fields, there exists an infinite number of possibilities to 
combine magnetic and electrostatic fields in order to increase one of the two 
abovementioned factors without any decrease of the other. Introductory theories 
treating such more general systems, have been developed among others by 
Cotte (1938) and WenpT (1942). It seems adequate, however, to select for a 
more close study some special groups, which are characterized both by simplicity 
and practical advantages. In the present paper the treatment is limited to field 
combinations possessing cylindrical symmetry. It is also aimed to be an in- 
troduction to a theory of a mass spectrometer for mass defect work which is 
being built at this Institute. An essential feature of this apparatus is the 
realization of almost ideal focusing conditions and the avoidance of fringing 
fields in the ordinary sense. 

In the mathematical treatment of the aberrations of such focusing systems, 
one has to choose between the eiconal method and the path equation method. 
Since in either case the cylinder symmetrical fields are unique to give explicit 
formulae, the author preferred the latter method. In more general cases where 
the aberration entities are obtained as integrals containing functions, which are 
defined only by differential equations, the eiconal method has evidently some 
advantages especially in the systematization of the aberrations. A thourough 
discussion of the latter ones is made by Marscuauy (1944) in a paper, which 
was brought to the author’s notice by Professor 8. Friiaar at the Birmingham 
conference in 1948. I wish to thank Professor FLUcan for his communication 
and Dr. Marscuatt for his kindness in sending me a copy of his paper, which 
has not been available here. Though aberration formulae for more general field 
forms are given, it seems as if just the eiconal method, which is used in that 
paper, has prevented the author to find the most interesting special focusing 
systems exhibited by fields of cylindrical symmetry. 
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2. The electric field 


The electrical as well as the magnetic field is assumed to have cylindrical | 


symmetry with respect to the z-axis and mirror symmetry with respect to the 
xy-plane. If the electrical field strengt’ is given within an annular region of 
the latter plane and is expressed by a series expansion in the vicinity of the 
circle r= 79 
,(r, 0) = E(r) = Ey (1+ Bio + fee? +--*) (2:1) 
where 
= (t— 1)/To 


one can compute the potential and, consequently, the field components G, (r, z) 
and €,(r,z) within an annular spatial region in the vicinity of the above- 
mentioned circle. Putting 


(1,2) = Dp (rn)? (2:2) 
we obtain ar 
0 S 
6, (r, arse (er, 2) =— >i g(r) 2” (2:3) 
v=0 
&, (7, =< v1, 2) =— D\y (r)2v22?71 (2:4) 
v=1 


and by means of the condition div © = 0 


1 0 f) 
FA AUS as 
Sf oy ’ 
an E At Gy (r) + (29 + 1) (29 + 2) p30 | Ze = (ee 
resale 
whence 
am 2 Pp gurg 
PO) = DIN > de” de) M0 (2:6) 


From Eqs. 2—6 we get the following explicit formulae 


Ss 


(r, 2) =— Eyrolo + $810" + $ 620° + tBs0* —3(1 + Bi) C+ 
+3 (1—B,—2 Ba) @ 0°—3 (1—fy + Bo + 3 Bs) 07 6? + sk (1—B, + 4 Bo + 6 By) C4] (2:7) 
G,(r, rani sls ase or ‘< + Bg Q* + 
3 (1 = Pi— 28s) Co — {1 > fa > Ba BBs) eC (ae 
(7, z= By [— (1 4981) o HB, = 2B oe (I By 4 Be Sie cee 
+ 6(1— Bi + 4 Bo + 68s) 0°] (2:9) 
where €=2z/r9. The potential is put equal to zero for T= 7h Zi =O, 
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For the purpose of realizing such a field we need also an explicit formula 
for the shape of the equipotential surfaces. Putting 


O(r,z) =— Lyra (2:10) 


_ where o is a constant, we may calculate 9 as a function of ¢?. By means of 
successive approximations we obtain 


e=o— 3 fo + (Shi — 5 f.)o° +43[1 + B.—(1 + Bi—2f,)o] 0%. (2:11) 
The curvature of a radial section of the equipotential in o = o,, ¢ = 0 then is 


lth (1 +t fi 2 fs) o, (2: 12) 


and the mean curvature in the same point 


Fy — fr — B- 2B) oy (2: 13) 


We may notice the special case 
2fs=1+ 8 (2:14) 


where all equipotentials have the same curvature. 


3. The magnetic field 


An analogous calculation of the magnetic field may be started from a series 
expansion of the axial field component in the xy-plane 


H.(r, 0) =H (r) = Ho (L + 01.6. 403.074 030% +%-). (3:1) 


Even in this case the potential of the field 


P(r, 2) = Dy pr (rer (3: 2) 
v=0 
and the field components 
D- (7,2) =— a-¥ (7, 2) =— Dy y(n) 274} (3:3) 
v=0 
0 < 
8. (7,2) = 5 (r.2) =— Dior) (Qv + 2” (3:4) 


are calculated for an annular region in the vicinity of r=79, 7 =0. In this 
case the condition div ) = 0 gives 
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Yo (r) =— H (r) (3:5) 


1 Ta ae ; 
vO =aroail- sata Yo (7) (3:6) 


and finally 
Y (r, 2) =— Horo (6 + ae6 + age? + age? —4 (a, + 2ay) C0? + 
+ 1 (a, —2a,—6a3)067] (3:7) 
Gy (r, 2) = Hy [aro + 2a,0f + 3ase76 + 4 (a, —2a,—6a,)C*] (3:8) 
Gz (r, 2)=Ho [1+ ay 0 + ag 0? + a3 GE? —3 (a + 2 ag) 2? +4 (ay—2 ag—6 a) 0 C7]. (3:9) 
The field components can also be derived from the vector potential 


A (r, z) = Hor [oe — 3 (1 — a) @® + § (38 —ay + 2 ag) 0? —75 (6—2 ay + a,—3 a) = 


— $0, 0° — a, 00° Page 0 Site 20, — Cac (3: 10) 


by means of the relations 


¢ A (r, 2). (3: 11) 
As in the electrical case the equations of the equipotentials are calculated 
from the scalar potential. We put 
Pt.) =— bathe (3: 12) 
and compute € as a function of @ by successive approximations 
€=06+ h(a, + 2a.) 07 — a, ce + (aj — ag) o 0”. (3: 13) 


This equation also determines the shape of the pole surfaces if iron of large 
permeability is used to produce the field. 


4. The equations of motion and their solution 


If a particle of mass m and electric charge e is moving in an electric field 
© and a magnetic field §, the laws of motion in cylindrical coordinates are 
expressed (in non-relativistic approximation) by 


m(r¥ — 7g") =—rpehz(r,z) + eG, (r, 2) (4:1) 
mi(rg + 2ry) =reHz(r, 2) —2eH,(r, z) (4 : 2) 
mi=rpeHr(r, 2) + eG, (7, 2). (4:3) 
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Eq. (4:2) can be integrated once with respect to the time coordinate ¢: 
po =~ rA(r,2) +0 (4:4) 


where A (7,2) is the vector potential of the magnetic field defined by Eq. (3:10). 
The integration constant C is determined from the initial conditions: the 
particle is assumed to proceed from the point r=79, z=0, gp =O with the 
velocity v, and with the inclination y from the tangent of the circle r= 7%, 
z= 0 in the same point. This gives 


C = 1901 COS y = 79 Uy. (4:5) 


Now ¢ is eliminated from Eqs. (1) and (3) and ¢ is substituted by » as an 
independent variable by means of Eq. (4), for example 


d 


dp 
a PAU SE el? . 
=¢ 75% ig7 yr Berle (4: 6) 


Thus the equations of motion are reduced to 


r’’ + B-1 pee +2 ) (ppt Ge Bees (4:7) 
Or Oz m m 
z’ + Bo (2 eS) of Rate OIG Bate © (4:8) 
0 Oz m m 
where 
e 1. WU 
Brea a0 (4:9) 


It is seen from Eqs. (5), (7) and (8) that the circle r= 79, z =0 is a possible 
trajectory if the velocity 1, has a value vj determined by 


%9 =) Bs 1. sally aa Bh ~ Eo (4: 10) 
and 
eA ieee 20. : 
By = etal: (4:11) 
or 
2 
5 = en Hy —e Ep. (4:12) 
, 


The latter equation expresses the balance between the centrifugal force and 
the magnetic and electrostatic forces and may be considered the zeroth ap- 
proximation of the solution of Eqs. (7) and (8). In order to obtain higher 
approximations both membra of the same equations are expanded in powers of 
0 =(r—7)/r and €=z/ro giving 
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0 + Apo Fo t+ Fe + Hg (4:13) 


t+ C=O + Gat Gs ter (4:14) 


where the first order terms in @ and € are assembled in the left membra, and 


all other terms in the right membra. The functions 7 and G may also be 
developed into series with respect to the velocity difference 0+vp9 = v1 — V% 
and the angle y. One should notice here that the functions are of the first 
order in 6 and the second order in y at least. 


5. Directional focusing 
a) The focusing angles 


The characteristic feature of the geometrical focusing is most easily explained 
by first putting v; =v», 7. e. assuming that a particle of velocity v9, determined 
by Eq. (4:12), emerges from the point r=79, z= 0 on the central line with 
an inclination y to the tangent of the central line. in the same point. If terms 


of higher order than the first in @ and € and y are neglected, we obtain after 


some calculation from Eqs. (4:7) and (4:8) 


oe’ +4e=0 (5:1) 
rere (5:2) 
where 
Ap ae 1) Ge 2 ee ae ee) (5:3) 
dg a Oe (5: 4) 
and 
# = ery Hy|mr (5 735 


is the ratio between the magnetic force and the centrifugal force. 
Eqs. (1) and (2) show that @ and ¢ are periodical functions of the angle p 
if A, and A, are positive. In that case the particle will pass through the 


cylindrical surface r=rg after an angle O” = an/ Var (n = 1, 2, (¥")"emm 
through the symmetry plane z = 0 after an angle 0” = xm Vae (m = 1, 2,...)s 
It follows that the so-called focusing angles 


@, = O© = x[1 + (x — 1) (x — 2) + ax + Bi (1 —x)]-# (5: 6) 
©, = © = a[x —1— a,x — B,(1—)]-* (5:7) 


only depend upon the inhomogeneities of the magnetic and electrostatic fields 
and the ratio x. The planes y= ®, and » = ®, are called focal planes for 
radial and axial focusing respectively. The simple relationship 


lees) lconglade CRUE 
DF nO pis me 


(5:8) 


shows that the focusing angles may have all values between w[1 + (x — 1)*]7! 
and co, A smaller value for one direction is always accompanied by a de- 
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focusing in the other direction. For double directional focusing, ®, = &, = @, 
the largest possible focusing angle ® = aV2 occurs at x =1, that is a pure 
magnetic field. 


b) Radial focusing 


The condition for radial focusing ®, = co or 2, = 0 gives according to Eq. (4) 


% = (1 + B)/(lL—a + f:) (5:9) 
or, by means of Eq. (5), 


€ Vp Hy|e Ey = x/ (% — 1) = (1 + By)/ ay. = (510% 


This result includes the following special cases some of which are well known 


a) x=1, a, =0, pure magnetic field, homogeneous and having the focusing 
angle ©, =z. 

6) x =0, pure electrostatic field with 6, = — 1, 2.e. the radial field between two 
cylindrical and concentric conductors. The focusing angle is ®r = 2/V2 ~ 127°. 
This focusing principle was proposed by Hucues and Rosansxy (1929) and has 
applications within different fields especially in precision mass spectroscopy. 

y) x = 2, simultaneous magnetic and electrostatic fields whose gradients are 
coupled by the condition 6; =2a,—1. The focusing angle is 6 =2/V2. The 
subcase a, = 0, 6, =— 1, that is homogeneous magnetic field and radial electric 
field was proposed by Bartky and Dempster (1929). It has found applications 
in the determination of e/m for electrons by SHaw (1933, 1938) and in a mass 
spectrometer constructed by Bonpy, JOHANNSEN and PoppER (1935). A char- 
acteristic property of all field combinations for which x == 2 is the first order 
velocity focusing, which will be discussed in Section 6. 

6) x ~ co, This case corresponds to an infinitely small focusing angle or to 
a rectilinear central line. The special case a, = 0, 6; =— 1 corresponds to the 
well known trochoidal motion. 


c) Axial focusing 


Since there seems to be no obvious applications of a pure axial focusing, 
we are content with an enumeration of the special cases corresponding to the 
cases a, 6, y for radial focusing. 


a) x=1 gives o, =— 1 and =a 
B) ~»=0 gives Bj =—3 and ©,=2/V2 


y) *=2 gives Bj =2a, +1 and @,=2/V2. 


d) Double directional focusing 


The condition for equally strong focusing in the radial as well as in axial 
direction gives according to Eqs. (5:6) and (5:7) 
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$,=0,=0=2| > |’ (5:11) 


and 
a,x + By (1—x) +4 (x— 2? =0. (5: 12) 


We notice the following important subcases. 


a) x=1, which gives a, =— 4 and 6 =2V2 ~ 255°. This case was applied 
in the large B-spectrometer constructed by K. Siegbahn and the present author 
(1946, 1949). The magnetic field strength decreases outwards in such a way 


that its first derivative is the same as for the function 1/Vr. 

B) x=0, which gives Bj =—2 and O=az. This case is the electrostatic 
counterpart to case a). An example of such a field is the electric field in a 
spherical condenser. Cf. PuRcELL (1938). 

y) x =2, whence 6j = 2a, and ®©=-z. As in the case of radial (and axial) 
focusing there exists an infinite number of field combinations, which also give 
velocity focusing of the first order. Among these cases we observe especially 


1) a, =/,=0, which means a homogeneous magnetic field combined with a 


radially constant electric field (not a cylinder field), 2) a, =—+4, Bj =—1, 
which is an electric cylinder field combined with a magnetic field of the type 
1/Vr. 


Since one of the quantities a, and 6, may be fixed arbitrarily there is one 
degree of freedom left for reducing the aberrations. This possibility will be 
discussed in Section 8. 


6. Velocity focusing 


We now make the more general assumption that the initial velocity 2, 
differs by 6:vp =v, — vp from the velocity vg determined by Eq. (4:12). Series 
expansion of the functions Fy and Gy» in powers of y now gives 


Ser CO toe eo) 


foment ag ay (1 + 6) 


2 
JF? G=0 (6:1) 


where terms of higher order than three are neglected. If only first order terms 
in ge, € and y are considered, we obtain the following differential equations 


é (6: 2) 
Ci dat =e 0, (6: 3) 


In general the coefficients 7, and A, are functions of 6 and y, but if terms of 
higher order than the first in y and 6 are neglected, we have the formulae 


Ay = 1, + 6[—8 + ll x — 42% — xa, — 2(1— x) By] (6: 4) 


Ae = de + 6(2— 2 He eye QL) By (6: 5) 
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where 4, and A; are given by Eqs. (5:3) and (5:4). The complete solutions 
of Eqs. (6:2) and (6:3) under the above-mentioned initial conditions are 


Y é = 62—x+6 = 
siete now 0 fy gears ; 
0 Vi cos 8 sin y VA, + i, + oF (1 — cos e Vir) (6: 6) 
aoe ; = 
¢ =— = sin 0 sin py VA, ar 
Vin . set 


where @ is the angle between the osculating plane of the orbit and the sym- 
metry plane z—0 in the starting point g=0. In the focal plane for radial 


focusing y =2/V42, 0 has the value 
0, = [(2 — x) + (2% — 3) 6] 26//, + 


+ F[8—1lx + 40? + way + 2(1—x) Bil y 4;# cos 8 (6 : 8) 


where terms of higher order than the second in 6 are neglected. This equation 
shows that the focusing system has a dispersion proportional to (2 — x)/A,. 
In the important special case x = 2 the linear dispersion disappears and Eq. 
(6:8) simplifies to 


0, = a(1 + a — B,) y bA-? cos 6 + 267/34, (6:9) 


This result can also be expressed in the following way. If the electric force 
just balances the centrifugal force, then the dispersive effect of the electric 
field eliminates the dispersion of the magnetic field. Formula (6:9) was derived 
in the special case aj = 0, By =—1, 4, =2, 6=0 by Barrxy and DempstTER 
(1929). 


7. The geometrical aberrations for a point object 
We return to the assumption v; =v) and compute the functions Fy and Gy 


defined by Eqs. (4:13) and (4:14) by a further expansion of both membra. 
of Eqs. (4:7) and (4:8) and obtain 


Fo == (= 2) 0" 10 ne [aac ie) 
where 
fia =(%—2) (2 x? 2 + 3) 4 AG 4) oy —2 (%—1) (x—2) By —x% ag—(1—x) By (7: 2) 
faa = 4 (% —1) + 5(e— ay + 31%) + Hag + (1%) Bo (7:3) 
and 
Gy =—(% — 2) 0' C + 91206 (ioe) 
where 


Bes =e — 1) (Die 3), +4 (3-2) ay (41) (24D) By 2% ag + 2(1—x) Ba. (7:5) 
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The differential equations | 
o'+A0=Ff, + Fs (7:6)4 


CY +426 = Gy + Ge (72 
where now 
Fy=3(%—2)y", Gy =0, (7:8) 


are solved by a method of successive approximations. In the first approxima- 
tion where the right membra of Eqs. (7:6) and (7:7) are put equal to zero 
the solutions are . 


(= cos 6 sin vy VA, (7:9) 
(0) KS: ; ~ 
Cs Vi. sin @ sin Viz (7: 10) 


where the initial conditions 
oh) 0; FO <0; 0 = pices 0-0 O== y sin OMor D = 0 (Tg 


are adopted. Next approximation is obtained by substituting the unknown 
functions @ and ¢ in F, and G, by the first order solutions (7:9) and (7:10), 
From the inhomogeneous differential equations so obtained we determine par- 
ticular solutions, which satisfy the following initial conditions 


a) aa) o) CY’ —0 for p= 0 (7:12) 
whence the complete solutions of the second order are 
= 0 ae 0, C= CO 46M, (7:13) 


By means of elementary calculations we find 


ot) i | cos? 6 + (. — 2+ fs) sin? | $y? (1 — cos pVA;) + 


dz 
1 fia 1 ,,2 2 4 fo 
anaes mt do 142 cos? 6 (cos 2y Var — cos pVA;) + 
fe eee Var) :14 
Aaa ee sin? 6 (cos 2 p VA, — cos 9 VA; (7: 14) 
: I = = 
oo Vie VEDA (« 2 pa) 14% sin 6 cos 6 [cos g (VA, — VAz) — 
1 g 
—vcos » VA,| + - 94 ee 
Vira = AG VA, ;) 


1? sin 6 cos O[cos p (VAr + Vaz) — cos oVaal. Ae 15) 
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It should be noticed that these formulae are applicable also in the limiting 
cases J, > 0 and A; > 0 as well as in the case of resonance J, > 4A). 

Since we are especially interested in radial focusing we compute the line 
width in the focal plane for radial focusing a putting g = 2/ VA, in Eqs. 
(7:14) and (7:15) 


0. = = | cos? 6 + (. —2+ fi) sin? | y+ sj (. 2 + A) v* cos 6 + 
r4 (0 Ip 


‘ty sin? 6[cos 2 VAz/A» + 1] (7: 16) 


2, 2 F —— 
Ce f (: 2 i) y* sin 6 cos 6 cos tt VAzl Ap. CELI 


In the case of two-directional focusing, these expressions are simplified into 


4 = 2 4 ; 
0 = (Sahn a \? cos! 0 + (2 + Zia), 9 sin? 6 = (7: 18) 


4 2 2g : 
at) ee 12 2 Pp 
SI oes = (» 2 4 )y sin 6 cos 6. (ee aR) 


It remains to give the explicit formulae in the important subcases, which 
were mentioned in Section 5. 


Aa) A= 1 aS = 


0) =— 2 (8 a, — 1) y” cos? 6 + 3 (8 a, — 3) y” sin? 6. (7: 20) 

6D = $ (8 ag — 3) y* sin 6 cos 6, (7:21) 
ty hy ee sg 258) 

of) =— 3 (2B, — 3) y? cos? 0 + $(B. — 3) y* sin? 0. (7229) 

CY = § (8, — 3) y? sin 4 cos 0, (7 : 23) 


Vy) 2=2, A=1-$,=24,: 
oP =— §(2 a, + 2a, — Bp.) y” cos?.6 + 2(1 + 4a, —2 fo) y? sin? 6. (7: 24) 
CY) = 4 (1 + 4a, — 2B.) y” sin 6 cos 6. (7: 25) 


8. Directional focusing of the second order 


a) Two-directional focusing 


Since high resolving power at a given energy or mass dispersion is a desir- 
able feature of every focusing system, we are mainly interested in the radial 
line width in the focal plane for radial focusing. Or, the quality of the radial 
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focusing is a factor determining both resolving power and transmission power 
of the system, while the axial focusing is important only to the transmission 
power. From that point of view we have to find all systems in which the 
radial line width o{) determined by Eq. (7:16), vanishes for all values of 6. 
According to Eqs. (7:2), (7:3) and (5:12) we then obtain for all values of x 
except 0 and 1 


6%? —8x+4+3 

=—14 Sia 

ay 4 (x ) 4x (x — 1) ( ) 
8x*—15%+8 

bey , 8:2 

By 3 + (x 2) 4(x—1)? ( ) 

Ay Bo. *#+2 x — 2 x 
= any 17 — : Sis 
lane xe Be BMA ae 


These equations show that second order focusing occurs at a suitable choice 
of the coefficients a,, 61, G2, Bz for all values of x except x = 1 (purely magnetic 
field) and x=O (purely electrostatic field). In the case of velocity focusing 
x = 2 we get 

a =, bp =, Be = 2a + 5. (8: 4) 


From Eq. (7:17) it follows that the axial line width ¢“ also vanishes in this 
case. Consequently every field combination fullfilling the conditions (8: 4) 
represents an almost ideal imaging system possessing both first order velocity 
focusing and second order two-directional focusing. 

b) Astigmatic imaging systems 


If we are not interested in axial focusing there might exist other fields, 
which provide only point-to-line representation but still second order radial 
focusing. An investigation of the case x = 2 shows that there are two groups 
of such cases, namely on one hand 


A=1—2a, and Bo = 20, + 2a (8:5) 
giving 
A= 4G, and A, =2—4ay 
and on the other hand a root of the equation 
cot ag=1—4”¢ (8 : 6) 
where ¢ = VAz/ Ap = ®,/@, together with the condition 


BP, = 2a, + 2ay. (8:7) 


In the extreme cases x =—1 and x= O there are no similar solutions. In a 
B-spectrometer constructed by Langer and Cook (1948) where the field has the 
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coefficients a,=0, a,—=— } the line width oe is zero only for paths in the 
‘symmetry plane while all other paths (0 ~ 0) give a line width 


of) = (2 — $x?) y? sin? 6. (8:8) 


This error, however, should be negligible as the axial aperture must be kept 
small because of the absence of axial focusing. 
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